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Oh! Abstract 

J|h ! Deviations from expectations have been claimed for solar, atmospheric and 

high energy prompt neutrinos from charm decay. This information, supplemented 



only by the existing very good upper limits for oscillations of the at accelerator 
energies, is used as input to a phenomenological three-flavour analysis of neutrino 
mixing. The solution found is unique and completely determines the mass eigen- 
states as well as the mixing matrix relating mass and flavour eigenstates. As- 
suming the mass eigenstates to follow the hierarchy m\ <C m 2 <C m 3 , their values 
are found to be m 1 < 1(T 2 eV, m 2 = (0.18 ± 0.06) eV, m 3 = (19.4 ± 0.7) eV. 
These masses are in agreement with the leptonic quadratic hierarchy of the see- 
saw model and large enough to render energy-independent any oscillation-induced 
phenomenon in solar neutrino physics observable on Earth. This possibility is not 
excluded by the present knowledge of solar neutrino physics. The mixing angles 
are determined to be 12 = 0.55 ±0.08, 6 13 = 0.38 ± 0.06 , 6 23 < 0.03. Small values 
of 623 are typical of any solution in which m 3 lies in the cosmological interesting 
region. The solution found is not in serious contradiction with any of the present 
limits to the existence of neutrino oscillations. The most relevant implications in 
particle physics, astrophysics and cosmology are discussed. 
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1 Introduction 



In recent years, several deviations from expectations have been reported in various 
fields of neutrino physics 

At high energy, the equality between the v e and energy spectra prescribed 
by e-fi universality for the semi-leptonic decays of charmed particles appears to 
be violated at the four sigma level. 

In atmospheric neutrinos, the ratio between the and v e fluxes is found to 
be substantially less than the predicted value of 2. 

Even if quantitative conclusions depend on the somewhat different predic- 
tions of the various solar models, measurements of solar neutrino fluxes on Earth 
have been known for some time to fall short of the corresponding theoretical 
expectations. 

These positive results have been taken as indications of the existence of non- 
zero neutrino masses and mixings. However, they have to be examined for con- 
sistency with each other and confronted with the many negative searches for 
oscillations and their corresponding upper limits. 

No satisfactory purely phenomenological over-all interpretation has been ob- 
tained so far. 

In this paper we describe a comprehensive three-flavour analysis which results 
in an unique, overconstrained and complete description of neutrino mixing. Neg- 
ative results being usually available only at the 1.64 sigma level (90% CL), we 
have used as input data only the three positive results above and the two most 
stringent upper limits obtained in the searches for v e -v^ and v^-v T oscillations at 
high energy accelerators. We also discuss the compatibility of this solution with 
the other existing limits and its most relevant implications in particle physics and 
astrophysics. 

Preliminary accounts of this work have been presented elsewhere ]2j, |3| . 
CP-invariance is assumed throughout the paper. Also, unless otherwise stated, 
the term neutrino is used to indicate both neutrino and antineutrino. 

2 Three-flavour formalism 

In the complete three-flavour approach, the weak eigenstates \u a ) = u e , u^, v T and 
the mass eigenstates = ^1,^2,^3 are related by an unitary transformation 
matrix U, in terms of which the probability of an initial neutrino v a of energy E 
being equal to another neutrino vp at a distance L, can be written as 

P a p = <W - 4 ^ U ai y pi U aj U Pj sin 2 (A 4J /2) (1) 
with Ajj = SmfjL/2E, where Sm^ = m? — m 2 , m 8 = m{i>i). 
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The [/-matrix can be parametrized as 

(C12C13 •S 12 Ci3 s 13 \ 

-S12C23 - C12S23S13 C12C23 - S12S23S13 S23C13 (2) 
S12S23 - C12C23S13 -C12S23 - S12C23S13 C23C13 / 

with Cij = cos 9ij and Sy = sin 6^, where #12, #13 and 6*23 are three independent 
real angles lying in the first quadrant. 

Of the three £m 2 's appearing in eq. ([!]), only two are independent. The 
complete solution of the problem consists therefore in determining five unknowns: 
two <5m 2 's and the three %'s. 

For rrii ^ m 2 ^ m 3, only two <5m^-'s characterize the oscillatory behaviour of 
eq. ([I]): 5m\ 2 and 5m\% — 5m 2 ^ 3> <5m 2 2 , corresponding, respectively, to a "slow" 
and a "fast" oscillation. Thus, for any given experimental situation, depending 
on the range of L/E under study, either only the fast or both the fast and slow 
oscillations may be occurring. 

This allows to define the two ranges l/5m 2 3 (eV~ 2 ) <C L/E(m/MeV) <C 
l/5m 2 2 (eV~ 2 ) ("short-baseline") and L/£(m/MeV) > l/<5m 2 2 (eV~ 2 ) ("long- 
baseline" ) in which the average transition probabilities P^p and P£p are calculated 
from eq. (0) for sin 2 (A 12 /2) = 0, (sin 2 (Ai 3 /2)) = 0.5 and (sin 2 (A 12 /2)) = 0.5, 
(sin 2 (A 13 /2)) = 0.5, respectively. 



3 Input data 

3.1 Accelerator neutrinos 

In prompt neutrinos from charm decay, the equality between the u e and spec- 
tra prescribed by e-fi universality appears to be violated as the neutrino flux 
asymmetry 

A = (u^ flux — v e flux) / {yp flux + v e flux) (3) 

is experimentally determined to be A = 0.21 ± 0.05 [|J. 

For values of L/E smaller than about 0.1 m/MeV, and consequently in the 
domain of the prompt neutrino data, P efl and P^ T are experimentally known to 
obey the 90% CL upper limits [H 

P eix < 1.5 x 10" 3 , P Mr < 2 x 10" 3 . (4) 

Thus, ascribing the non-vanishing value of A to a depletion of the v e flux 
due to oscillations, the analysis of the L/E dependence of the data leads to the 
existence of v e -v T oscillations with the parameters ||, |(| 

sin 2 (2a) = 0.48 ±0.10 ±0.05 (5) 

5m 2 = (377 ±27 ±7) eV 2 . (6) 
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Table 1: Atmospheric neutrino experimental results. 

EXPERIMENT r = (/Ve) obs /(/Ve) calc 

Nusex 1.04 ±0.32 

Frejus (contained events) 0.87 ± 0.19 

IMB-3 0.54 ±0.13 

Kamiokande 0.60 ± 0.08 

Soudan 2 0.64 ±0.19 



In terms of three-flavour notations, eq. (0) ensures this to be the fast oscillation 
(sin 2 (Ai 2 /2) = 0), and the limits of eq. @ imply that 6*23 must be very small. It 
then follows that eqs. (|) and (|6|) can be rewritten as 

sin 2 (2#i 3 ) = 0.48 ±0.12 (7) 

5m 2 13 ~ 5m 2 3 = (377 ± 29) eV 2 . (8) 

3.2 Atmospheric neutrinos 

Neutrinos are also produced in the interactions of the primary component of 
cosmic rays in the Earth's atmosphere and in the subsequent decays of the pro- 
duced secondaries. As almost all these decays involve a muon, which in turn also 
decays, the ratio R between the and v e fluxes can be safely predicted to be 
approximately 2. 

Several experiments have addressed this question f/J. They all consistently 
find values of R which are smaller than 2 by a factor of r, approximately equal 
to 0.6. In fact, for average neutrino energies below 1 GeV, the weighted average 



of all available data 0, [|, [10], [TTJ, [T^] listed in Table 1 gives 

r = 0.63 ± 0.06 with xVd-o.f. = 0.96. (9) 
Contrary to the low-energy behaviour, in the multi-GeV energy range r is 



observed to be a function of the zenith- angle ||13||. This is very suggestive of an 
L/E dependence of r and thus of an oscillation- induced phenomenon. The 5m 2 
characterizing this oscillation lies in the 90% CL interval 

5 x 10~ 3 < 5m\ 2 < 8 x 10" 2 eV 2 . (10) 

Eq. (|1(J) ensures that in the sub-GeV energy range all detected observables 
have their average long-baseline values. Thus, in terms of the three-flavour for- 
malism, eq. (H) becomes 

{P^ + PPe,){P' l Pe, + Pee)'' = 0-63 ± 0.06 (11) 

where p = 0.47±0.02 is the expected v e jv il flux ratio in the absence of oscillations 

MM 
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Table 2: Solar neutrino experimental results and theoretical predictions (lcr). 
Statistical and systematic errors have been added in quadrature. For the Ga 71 
entry, the quoted experimental value is the weighted average of the SAGE (69 ± 
" SNU 



10t? SNU H) and Gallex (77.1 ± 8.51^ SNU ||) results. The error on this 
value is almost certainly underestimated as systematic errors have been treated 
as uncorrelated. 



DETECTOR 



Ga Yi 
(SNU) 



CI 



Kamiokande 



(SNU) (10 6 cm-V 



Experimental result 2.55±0.25 2.75±0.45 73.7±7.6 
Princeton/Yale model 8±1 5.7±0.8 132±6 

Saclay model 6.4±1.4 4.4±1.1 123±7 



3.3 Solar neutrinos 

Solar neutrinos have been known for some time to deviate from expectations. 



The experimental results ]T6| [L7|, [18], 19 , coming from three different types of 



detectors, as well as the predictions of the Princeton- Yale |^0| and Saclay [^] solar 
models are summarized in table 2. With the only exception of Ga 71 , experimental 
results are dominated by systematic errors. Furthemore, it can be seen that 
theoretical uncertanties almost always exceed experimental errors. Thus, in any 
comparison between models and experiments, theoretical errors must necessarily 
be taken into account. As they are highly correlated, we have limited our choice of 
theoretical models to those published and for which the complete error correlation 



matrix is available [22 



With neutrino mass differences as large as those of eqs. (§) and fllQ|) , all solar 
neutrino oscillation-induced phenomena observable on Earth are bound to be 
energy-independent. In particular, in all experiments the v e flux must be reduced 
by the same factor F, independently of the detection threshold. 

Although marginally in some cases, the data are consistent with this expec- 
tation. Fig. 1 shows the two % 2 's obtained by comparing the experimental data 
with the theoretical predictions of the Princeton- Yale (PY) and Saclay (S) solar 
models as a function of F. In the analytical formulation of these x 2 ' s > theoretical 
and experimental errors have been added in quadrature and correlations among 
theoretical errors have been fully taken into account by means of the appropri- 
ate correlation matrices. The Xpy an d Xs minima turn out to be Xpy min = 9-48 
(CL = 0.9%, equivalent to 2.6 sigma) and xl ■ =5.10 (CL = 7.8%, equivalent 
to 1.8 sigma). The confidence level of the first fit is admittedly rather low, but 
not unacceptable. In view also of the goodness of the second fit and of the only 
approximate knowledge of both theoretical and experimental errors (particularly 
with regard to their Gaussian-like behaviour), no conclusion about a compelling 
energy-dependence of the F-factor can be drawn at this stage. 
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Figure 1: The functions x 2 = X 2 {F) obtained by comparing the experimental 
data with the predictions of the Princeton- Yale (dashed line) and Saclay (full line) 
solar models (see Table 2). F is the energy- independent factor by which all v e 
fluxes predicted by the same model are reduced. For each model the correlations 
among errors on the expected rates for the various detectors have been taken into 
account through the use of the appropriate error correlation matrix. All errors 
have been assumed to be Gaussian-distributed. 



6 



Following the procedure normally used by the Particle Data Group ||, the 
errors on the F-factors have been multiplied by the appropriate scale-factors. 
Then, the two F-factors turn out to be F PY = 0.46 ± 0.13 and F$ = 0.57 ± 0.13, 
implying the existence of a solar neutrino deficit. 

In order to quantify this deficit for our further analysis, we have calculated 
a new F-factor taking the arithmetic average of Fpy and F$- We have taken 
the error on this average (±0.08) to be a good estimator of a systematic error 
reflecting the fact that the two models of the same Sun yield different predictions. 
We have added this error in quadrature to the typical statistical error (±0.13), 
arriving at the result 

F = 0.51 ±0.15. (12) 

Using the formalism of the three-flavour analysis, eq. ( [12] ) trivially translates 
into 

P e L e = 0.51 ±0.15. (13) 



4 Mass eigenstates and mixing matrix 

Taking m-i -C rri2, eqs. (|8[) and (|T0| ) yield 

m!<10" 2 eV, m 2 = (0.18 ±0.06) eV , m 3 = (19.4 ± 0.7) eV (14) 

These values are consistent with the predictions of the see-saw model |23| . 
Using a leptonic quadratic hierarchy and the above value of m 3 one has mi = 
1.6 x 10~ 6 eV and m 2 = 6.8 x 10" 2 eV. 

From the three eqs. (0), ([TT|) and (|13"D, within the constraints (J|), the three 
angles 0^ are uniquely determined to be 

012 = 0.55 ±0.08 , 0i3 = 0.38 ± 0.06 , 23 < 0.03 . (15) 

with xLn = 0-093. 

As in any solution with m 3 larger than a few eV, the smallness of 023 is dictated 
by the upper limits and this implies an over-determination of the angles 0i2 
and 013. The good consistency of the input data is illustrated in fig. 2, which 
shows the three relations between i2 and 0i3 obtained from eqs. (|7]), ([TT]) and 
( |T3"D for 23 = 0. It can be seen, for instance, that the value of 0i3 is determined 
not just by the result of eq. (0) but also by the system formed by the other two. 
The ^/-matrix is then 

/ 0.79 ±0.05 0.49 ±0.06 0.37 ± 0.05 \ 
U= -0.52 ±0.06 0.85 ±0.05 < 0.03 I (16) 

V -0.31 ±0.04 -0.20 ±0.05 0.93 ± 0.03 / 

The knowledge of the angles 9^ allows also to calculate the following short- 
and long-baseline average transition probabilities P^g and Pfg 

Pg = 0.77 ± 0.06 P^ < 3 x 10~ 4 Pg = 0.23 ± 0.06 

> 0.998 P£ < 1.5 x 10- 3 

PL = 0.77 ±0.06 
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Figure 2: The three relations (±lcx) between 9 12 and #13 obtained from eqs. ([7|) 
(full lines), (0) (dashed lines) and (0) (dotted lines) for # 23 = 0. 
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P e L e = 0.47 ± 0.06 P e L = 0.35 ±0.05 P e L T = 0.19 ± 0.05 
p£ = 0.60 ± 0.06 P^ = 0.05 ± 0.03 

P T L r = 0.76 ± 0.06 

All upper limits are at the 90% CL. 



5 Discussion 

The results for the mixing angles (eq. (^)) are that 812 is fairly sizeable, 823 
is small and #13 lies somewhere in between the two. Small values of 823 are 
typical of any solution in which the mass of the heaviest eigenstate (eq. (0)) is 
of cosmological relevance j23 . 



The main implications of this scenario are discussed below. 
5.1 Accelerators 

Short-baseline experiments like Chorus and Nomad |25| are expected to see a 
u T signal. If not from z/^-z/ r [P^ T is consistent with zero), at least from v e -v T 
transitions (owing to the about 1% v e component in the beam and the sizeable 
P s ). 

Long-baseline experiments |7|, |26|] have a better chance to detect v tr v T oscil- 
lations, but the largest effect is anticipated in the v e - v ^ channel. 

Experiments at low energy accelerators like KARMEN [ETJ and LSND 



ought to be able to detect Ve-v^ transitions as, owing to the onset of the slow 
oscillation, P eil in their accepted range of L/E is typically a fraction of a percent. 
Both experiments should be in the position of investigating the conservation of 
the v e flux, a study which so far has resulted in the KARMEN not yet statistically 
significant Pf e lower limit [ p^j . 



5.2 Beta decay and reactors 

For Majorana neutrinos, under very reasonable assumptions, the neutrino-less 



double beta decay amplitude is proportional to |f29 



m W = J2 U ei m i ( 17 ) 

i 

The value of mpp calculated from eqs. (|TJ]) and (|16|) is 2.6 eV, while its present 
upper limit is only about several eV [[J, so that experiments will soon be in the 
position of providing definite conclusions on this very important issue. 

In about 13% of all /3-decays a heavy neutrino of mass ?»3 (eq.(|T4"D) is expected 
to be produced. Two-state analyses of the spectra end-points |30| might be 
able to investigate this question. However, at the present time, the presence of 
ununderstood phenomena in the high energy region precludes any firm conclusions 
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All /3-sources are expected to yield lower-than-canonically-calculated v e fluxes 
(both P^ e and P^ e are smaller than 1 and the wave-length of the fast oscillation is 



about 5 mm/MeV) but the only result available so far [32] is not accurate enough 
to really test this possibility. 

This effect should also be present in experiments at nuclear reactors. After 
some initial tantalizing results indicating a large depletion of the v e flux at the two 
sigma level [33. |31[], the more recent experiments have failed to substantiate any 
deviation from expectations. In fact, their claims represent the most stringent 
limit to oscillations of the v e and the most serious challenge to the result of eq. 

(§• 

For mass differences as large as that of eq. (P), owing to the very short wave- 
length of the oscillation, these experiments consist in comparing the antineutrino 
flux expected from the reactor with that experimentally measured at some (short- 
baseline) distance away. In spite of the many such experiments carried out in the 
last decade [36], [37] [38], |39| , no significant progress has been made in the ac- 

39| (the two best 



curacy of the results. The comparison of ref. |3j| with ref. 
documented and, respectively, the oldest and newest papers) shows that in both 
cases the determination of the ratio between results and expectations is limited 
by a systematic error of about 6%. This error is largely dominated by a common 
uncertainty on the knowledge of the reactor antineutrino energy spectrum. This 
is due to the fact that the bulk of the data on thermal neutron fission induced /3- 
spectra [40, 41, |42| all reactor antineutrino energy spectra calculations are based 
on were obtained only once, more than a decade ago. It has been argued repeat- 
edly that such important measurements ought to be checked, that the procedure 
used in deriving neutrino spectra from /3-spectra ought to be better understood, 
that the quoted systematic uncertainty is probably underestimated and that the 
over-all normalization error in the reactor antineutrino energy spectrum is more 
likely to be around 10% [43], [44], f4~5|. In any case, even taking the quoted 6% 



at face value and assuming it to be Gaussian-distributed, the reactor result dif- 
fers from the value of Pf e resulting from eq. ([T^) by 2.7 sigma, indicating some 
discrepancy but not a really serious incompatibility. 

Reactor experiments are potentially sensitive also to the existence of the mass 
difference of eq. flit)]). In this case the wavelength of the oscillation is certainly 
long enough to allow a comparison between two measurements of the antineutrino 
flux in two detectors at different distances from the reactor, thus reducing the 
systematic effects discussed above. Because of the low event rate, experiments so 
far have been barely able to investigate the L/E region of interest. However, the 
claimed limits |35|, |37|, |39[ exclude in part or almost completely the region 

even if in a different context, the re- analysis 

(0.085 ± 0.010) eV 



of eq. (|T0|). On the other hand 
of the reactor data of ref. 
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has shown that a mass m 
is quite consistent with all existing evidence. Reactor experiments with an even 
slightly larger L than it has been possible so far are clearly needed to clarify this 
important point. Long-base experiments will also have the benefit of a projected 
larger depletion of the v e flux [P^ e < P^ e ). 
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In conclusion, although reactor experiments do indicate the existence of some 
potential problems at the very edges of the explored region in the sin 2 (2a), 5m 2 
plane, the large uncertainties involved preclude any firm conclusion at the present 
time. 



5.3 Astrophysics and cosmology 

Neutrinos with masses of the order of ten eV can be the dark matter particles 
needed to explain the high rotation velocities in the outer parts of spiral galaxies 



47], 58] and the observed X-ray emission from hot diffused gas in elliptical galaxies 



Heavy neutrinos lead also to some important cosmological consequences. With 
masses as large as those of eq. fll4]), neutrinos contribute far more than baryons 
to the total matter density of the universe Q = po/p c where p c is the critical 
density (the subscript indicates present-day values). 

In terms of the Hubble parameter H (H = 100 h kms -1 Mpc -1 ) the baryon 
density is known to be || 

Q b hl = (0.015 ± 0.005) (18) 
while the neutrino density fl u obtained from eq. fll4|) is 

VL u hl = 1.075 x 10~ 2 x^m 4 = (0.211 ± 0.008). (19) 

i 

Thus, for a universe whose only massive stable components are baryons and 
neutrinos, the total density is given by 

Q h 2 = ( n b + tt u )h 2 = (0.226 ± 0.009) (20) 



The relation resulting from eq. fl20|) is shown in the D,q, Ho plane of Fig. 3 by 
the full-line curves. They define the allowed region (la) for a universe composed 
only of baryons (visible matter) and neutrinos (hot dark matter, HDM). The 
allowed region for a possible third massive component (cold dark matter, CDM) 
lies on the right-hand side of these curves. The same figure shows also the relations 
between Qq an d Hq obtained from the Einstein's equations for the two limiting 
values of the age of the universe T = 13, 16 Gyr and for two different values of 
the cosmological constant A (A = 3Hq\ /c 2 ) A = 0, 1. 

To date, the value of H Q is not precisely known, the various experimental 
determinations ranging between the values ho — 0.5 and ho — 0.8 |50|| . 

Depending on the value of H Q , different situations can occur. 

1. For ho — 0.5 the choice A = is consistent with the accepted bounds 
on the age of the universe. In this case eq. (|20|) requires Qo ?t 0.85 and 
for Qq — 1 neutrinos turn out to be the most massive component of the 



11 




Figure 3: Relations between the total matter density of the universe f2 (in 
units of critical density) and the normalized Hubble parameter ho (H = 
100/i kni s~ 1 Mpc _1 ). The full lines are for a sum of the neutrino masses of 
(19.6 ± 0.7) eV. The other curves represent the function h = h (il ) for the two 
limiting values of the age of the universe T = 13 Gyr and T = 16 Gyr and for 
values of the scaled cosmological constant Ao (Ao = Ac 2 /3Hq) Ao = (dashed 
lines) and Aq = 1 (dotted-dashed lines). 
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universe. For only baryons and neutrinos contributing to the total density 
of a critical universe, the result of eq. (|20"D implies 



H = (47 .4 ± 0.9) km s" 1 Mpc -1 

(21) 

T = (13.8 ±0.3) Gyr. 

2. For h ~ 0.65, eq. (PH|) requires A > 0. For A = 1, ~ 0.5 and the pres- 
ence of cold dark matter could be accomodated with a ratio flcDM /^hdm 
(flcDM /&v) which, depending on the value of flo, could be as large as ~ 0.5. 

3. For ho ~ 0.8 and Ao = 1, — 0.4 and the mass of the universe is again 
due solely to baryons and neutrinos. Non-vanishing values of flcDM /^hdm 
are possible only for Ao > 1. 

Thus, for any reasonable values of H and Ao, neutrinos with masses as large 



as those of eq. flTil) are the dominant component of the mass of the universe. As a 



consequence, their existence contradicts all cold-dark-matter-dominated (CHDM) 



cosmological models ||51|| . To have the same success that these models have in 
explaining the present data on the power spectrum of density fluctuations, hot- 
dark-matter-dominated models will have to rely on some new hypotheses, such 
as, for instance, seeds of non-Gaussian origin K 



6 Conclusions 

This analysis shows that, at least taking results at face value, the evidence for 
the existence of neutrino oscillations is really substantial. Although none of the 
existing positive signals is admittedly beyond some criticism, their ensemble is 
quite compelling. 

The signature of neutrino oscillations lies in the sin 2 factor of eq. ([I]) and any 
convincing proof of their existence must rely on some experimental observation of 
an L/E dependence. This evidence, first obtained from the analysis of the L/E 
modulation of the v^-v e asymmetry in prompt neutrinos 0, [J, has been further 
strengthened recently by the observation of the dependence on the zenith-angle 
(effectively also L/E) of the ratio between the and v e fluxes in atmospheric 
neutrinos of multi-GeV energy . 

All the available experimental information from accelerator, atmospheric and 
solar neutrinos is accounted for in the framework of a three-flavour neutrino os- 
cillation analysis. The solution found is unique and not significantly contradicted 
by any existing result, all conflicting evidence being below the three sigma level. 

Neutrino physics is presently at a crossroads. If is indeed large enough to 
be of cosmological relevance, the limits of eq. (|j) imply that 623 is very small. 
The relations between # 12 and #13 are then those of fig. 2 and all indications are 
that #i3 is sizeable. Under this condition, high energy accelerator experiments 
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such as Chorus and Nomad should not fail to see a u T signal. If they do, this 
implies that at least some of the input data to our analysis are wrong and that 
either #13 or 777,3 or both are smaller than the results we have arrived at. In the 
first case, the atmospheric neutrino result (see fig. 2) implies the relative large 
value for the F-factor F = 0.70 ± 0.07, thus disfavouring solar models such as 
that of Princeton- Yale [20] which (see fig. 1) require F H 0.5. In the other two 
cases 7773 must be smaller than a few eV, thus excluding neutrinos as the main 
constituents of the dark matter of the universe and leaving them out of a primary 
role in astrophysics. 
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